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NUMERICAL  SIMULATION  OF  HIGH  ENERGY  LIGHT 
ION  BEAM  FLOW  THROUGH  A  FOILLESS  DIODE 


I.  INTRODUCTION 

Interest  remains  high  in  the  quest  for  a  compact,  high  energy,  intense  light-ion  beam  accelerator. 
Some  of  the  areas  of  application  for  such  an  accelerator  include  directed  energy  systems,  inertial 
confinement  fusion  drivers  and  magnetic  fusion  fuel  injectors.  The  desired  operating  parameters  for 
such  an  accelerator  are  an  ion  energy  of  from  0.1  to  1  GeV,  a  beam  intensity  of  from  10*^  to  10‘* 
ions/cm^  per  pulse,  and  a  pulse  length  of  from  less  than  one  nanosecond  to  a  millisecond.  To  be 
termed  "compact,”  the  device  should  not  measure  more  than  10  to  30  meters  in  any  dimension. 

The  ability  to  generate  intense  relativistic  electron  beams  (IREB’s)  in  foilless  diodes*'^  has 
brought  the  concept  of  a  compact  ion  accelerator  very  close  to  a  reality  in  the  form  of  collective  ion 
accelerators  (CIA’s).  These  collective  accelerators  have  received  much  attention  over  the  years.^"* 
Their  operation  relies  upon  the  existence  of  localized,  collective  field  gradients  in  excess  of  I  MV/cm 
in  an  IREB.  In  comparison,  conventional  accelerators  provide  electric  fields  on  the  order  of  only  0.01 
MV/cm  which  is  unacceptable  for  a  "compact”  device.  Although  many  approaches  exist  for  achieving 
the  acceleration  of  positive  ions  via  the  collective  electric  held  produced  by  an  intense  relativistic  elec¬ 
tron  beam,  not  all  approaches  can  be  scaled-up  to  achieve  high  ion  energies.  Scaieabiiity  implies  a  sys¬ 
tematic  procedure  for  maintaining  synchronization  between  the  ion  being  accelerated  and  its  accelerat¬ 
ing  electrostatic  potential  well.  This  requirement  of  scaieabiiity  may  be  generalized  to  any  type  of 
accelerator. 

The  "space  charge  wave  accelerator"  concept  has  been  proposed  at  NRL  to  achieve  such  scaleabil- 
ity.^  The  fundamenul  idea  behind  the  concept  is  the  sustained  matching  of  ion  beam  velocity  with  the 
phase  velocity  of  the  accelerating  IREB  space  charge  wave  along  the  entire  length  of  the  CIA  (collective 
ion  accelerator).  The  most  formidable  challenge  to  this  concept  is  the  matching  of  iitjected  ion  beam 
velocity  with  the  minimum  electron  space  charge  wave  velocity.  Technical  restraints  appear  to  place  a 
lower  limit  of  0.2-0.3  c  on  the  phase  velocity  of  the  fundamenul  mode  of  a  1  MeV  IREB  space  charge 
wave.  This  dicutes  an  injected  proton  energy  of  16-30  MeV.  If  that  synchronization  of  initial  veloci¬ 
ties  can  be  met,  the  geometry  of  the  accelerator  cavity  can  be  designed  to  keep  v  and  v  p^ue  in  step. 
The  specific  design  used  is  termed  the  Converging  Guide  Accelerator  or  CGA.*  The  phase  velocity  of 
the  wave  in  a  conducting  cylinder  depends  on  the  ratio  of  beam  to  wall  radii.  The  phase  velocity 
increases  inversely  with  this  ratio.  Thus,  the  wave  can  be  accelerated  by  accelerating  the  beam.  The 
beam  can  be  accelerated  by  passing  it  through  a  meul  tube  with  a  converging  radius. 

In  Fig.  1,  ions  are  injected  from  the  right  along  the  central  axis  of  the  electron  beam  which  is  also 
propagating  from  right  to  left.  Next,  a  space  charge  wave  is  excited  and  grows  around  the  ions.  The 
ions  loaded  into  the  wave  are  accelerated  by  it  as  it  too  accelerates  in  step.  The  electron  beam  is  then 
dumped  and  the  ions  extracted  from  the  far  left.  Two  meters  of  magnetic  field  at  a  strength  of  20  kG 
are  required  for  this  system.  Initially,  the  electron  beam  planned  for  use  with  this  system  will  have  an 
energy  of  0.5  MeV  and  a  current  of  1.4  kA  for  250  nsec.  Also,  the  electron  beam  generator  can  be 
operated  at  10  pulses  per  second.  The  trace  that  is  shown  in  Fig.  1  is  a  photograph  of  the  machine  vol- 
uge  versus  time.  The  high  voluge  machine  was  used  in  conjunction  with  a  resistive  load.  For  proof- 
of-principle  testing,  low-current,  30  MeV  proton  beams  could  be  injected  from  a  conventional  ion 
accelerator  such  as  the  cyclotron  at  NRL. 
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Such  a  low  ion  current  experiment  can  indeed  test  the  wave  accelerator  concept  but  there  remains 
another  question  to  be  answered  before  kiloamp-levei,  ii\jected  beam  devices  are  constructed.  What  are 
the  consequences  of  injecting  into  a  foilless  diode  a  30  MeV  proton  beam  of  comparable  amperage  to 
the  current  normally  carried  by  the  electrons  alone  in  the  diode.  How  will  the  structure  of  the  potential 
well  inside  the  hollow  electron  beam  be  effected?  Will  the  net  electron  current  flowing  through  the 
diode  increase  and.  if  so,  by  how  much?  These  are  the  matters  addressed  by  the  computational 
research  reported  herein.  The  findings  demonstrate  that  an  iitjected  30  MeV  proton  beam  carrying  10% 
of  the  normal  diode  electron  current  will  not  seriously  effect  the  potential  well  structure  in  or  the  elec¬ 
tron  current  of  the  foilless  diode.  A  five-fold  increase  in  the  ion  beam  current,  on  the  other  hand, 
destroys  the  well  and  doubles  the  electron  cunent.  Thus,  clearly  *safe*  and  clearly  "unsafe”  ion  current 
levels  have  been  established.  The  determination  of  an  exact  "critical”  ion  current  for  use  in  the  CGA 
device  is  not  important  at  this  time  and  has  been  left  for  future  research. 

In  the  section  that  follows,  the  actual  physical  device  to  be  simulated  is  described  in  detail.  Sec¬ 
tion  III  then  outlines  the  key  characteristics  of  the  DIODE2D  computer  code  which  was  used  to  con¬ 
duct  the  numerical  simulation.  The  numerical  results  are  presented  in  Section  IV  including  such  infor¬ 
mation  as  radial  profiles  of  beam  current  density  and  plots  of  the  elecuostatic  potential  well  structure  in 
the  diode.  In  order  to  give  added  emphasis  and  understanding  to  the  numerics.  Section  V  is  devoted  to 
a  simple  theoretical  analysis  of  the  physics  of  the  ion  beam  and  diode.  Good  agreement  is  found 
between  the  theoretical  predictions  and  the  numerical  findings.  Finally,  Section  VI  summarizes  the 
major  conclusions  that  can  be  drawn  from  this  work  and  suggests  material  for  future  research. 

II.  THE  FOILLESS  OIOOE 

The  physical  diode  designed  for  the  experimental  testing  of  this  new  CGA  concept  closely  resem¬ 
bles  that  shown  in  F''g.  2.  It  corresponds  roughly  to  the  "wave  growth”  section  of  the  accelerator  shown 
in  Fig.  1.  The  coaxial,  vacuum,  power  feed  line  tapers  inward  radially  as  it  traverses  the  axial  length  of 
the  near-conical  cathode.  The  tip  of  the  cathode  coincides  with  the  entrance  plane  of  a  long  (>1 
meter)  cylindrical  drift  tube.  The  outer  radius  of  the  cathode  is  about  1.5  cm.  A  hole  O.S  cm  in  radius 
is  bored  through  the  center  of  the  cathode  to  serve  as  the  ittjection  port  for  the  high  energy  ion  beam. 
The  coaxial,  conical  A—K  gap  has  a  length  of  about  five  centimeters.  The  drift  tube  is  given  an  inner 
radius  of  about  1.7  cm.  The  entire  apparatus  is  immersed  in  a  uniform  axiai  magnetic  field  of  about  10 
kilogauss. 

Ideally,  the  model  diode  to  be  numerically  simulated  would  correspond  exactly  to  that  depicted  in 
Fig.  2.  Unfortunately,  the  simulation  code,  DIOOE2D,  which  was  used  in  these  studies  was  incapable 
of  efficiently  handling  that  specific  geometry.  Particularly  troublesome  are  the  sloped  walls  of  the 
cathode  and  of  the  coaxial  power  feed  line  as  well  as  the  overall  meter-long  length  of  the  physical  sys¬ 
tem.  The  existing  code  was  capable  of  treating  only  conducting  surfaces  and  boundaries  which  are 
strictly  radial  or  strictly  coaxial  (i.e.,  they  must  have  a  rectangular  /I -Z cross-section).  Thus,  a  cylindr¬ 
ical  cathode  shank  of  some  finite  thickness  must  be  substituted  for  the  conical  one.  Furthermore,  a 
combination  of  limited  resolution  and  finite  computer  data  storage  capabilities  prohibits  the  accurate 
treatment  of  a  system  less  than  2  cm  in  radius  but  over  100  cm  long.  The  toul  axial  length  of  the 
model  system  will  be  shortened  instead  to  10  cm.  The  peculiarities  of  the  Poisson-solver  for  the  elec¬ 
tric  and  magnetic  potentials  in  that  system  also  required  that  the  tube  be  terminated  by  a  fixed- 
potential,  conducting  wail.*  The  final  approximation  to  the  original  system  took  the  form  shown  in  Fig. 
3. 


The  cathode  used  in  the  numerical  model  is  a  hollow  cylinder  1.2  cm  in  radius  with  a  shank  wall 
thickness  of  0.4  cm.  It  projects  axially  2.5  cm  into  the  lO.O  cm  long  drift-tube  (anode).  The  anode 
tube  inner  radius  is  at  1.7  cm.  This  leaves  a  radial  A~K  (anode-cathode)  gap  of  0.5  cm  which  very 
comfortably  exceeds  the  0.17  cm  gyroradius  of  a  0.5  .MeV  electron  in  a  10  kG  magnetic  field.  Analysis 
using  existing  theory'”''^  leads  to  estimates  of  foilless  diode  electron  beam  currents  of  8-10  kA  or  diode 
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power  levels  of  4-S  GW.  Along  the  central  axis  of  this  device  will  be  ipjected  a  0.7  cm  radius  beam  of 
30  VfeV  protons.  Three  separate  simulations  were  conducted  of  this  configuration.  The  first  was  a 
'benchmark”  run  to  determine  the  normal,  electrons-only  operating  characteristics  of  the  foilless  diode. 
Of  importance  is  not  oniy  the  gross  e-beam  current  but  also  the  beam  profile  and  the  structure  of  the 
electrostatic  potential  well  inside  the  beam.  The  final  two  simulations  were  conducted  with  the  high- 
energy  proton  beam  flowing  through  the  device.  Specific  ion  beam  currents  of  1  kA  and  S  kA  were 
chosen  for  these  two  test  cases. 

III.  NUMERICAL  IMPLEMENTATION 

A  modified,  electrostatic-magnetostatic  version  of  NRL's  D10DE2D  computer  code  was  used  to 
carry  out  the  simulations.  The  details  of  the  code  may  be  found  elsewhere.'^  It  is  sufficient  here  to 
note  that  the  code  is  2-1 /2-dimensional  in  that  it  explicitly  moves  particles  in  R  and  Z  (radially  and  axi¬ 
ally)  in  space  while  keeping  track  of  all  three  components  of  their  canonical  momentum.  The  treatment 
is  fully  relativistic  for  both  electrons  and  protons.  Unlike  a  full  electromagnetic  simulation  code, 
DIODE2D  is  designed  to  find  only  equilibrium  charge,  current,  and  field  configurations  in  the  device 
under  study. 

Particle  contributions  to  net  charges  and  currents  as  well  as  the  resulting  E  and  B  field  values  are 
stored  in  a  mesh  of  discrete  dau  points  covering  the  entire  diode.  The  spacing  between  data  points  was 
fixed  at  R  ~  O.OS  cm  radially  and  Z  ~  0.078  cm  axially.  This  allowed  for  34  data  cells  to  span  the  1.7 
cm  from  the  centerline  to  the  inner  radius  of  the  anode  tube  and  128  cells  from  the  back  of  the  cathode 
to  the  anode  endplate.  The  protruding  cathode  shank  measured  8  data  cells  thick  and  32  cells  long.  On 
its  front  face,  electron  emission  was  permitted  from  the  outer  7  of  the  8  possible  emission  points.  The 
eighth  cell  was  not  allowed  to  emit  in  order  to  obtain  additional  numerical  accuracy  in  the  field  values 
governing  electron  emission  from  the  inside  lip  of  the  cathode  shank.  Emission  was  permitted  all  along 
that  inner,  as  well  as  the  outer  shank  surface  except  for  the  rearmost  6  cells. 

The  magnitude  of  the  electron  emission  from  a  given  cathode  surface  is  determined  by  the 
amount  of  space  charge  necessary  to  zero  out  the  perpendicular  component  of  the  electric  field  at  the 
emission  point.  In  order  to  simulate  the  injection  of  the  30  MeV  proton  beam  of  predetermined 
amperage,  the  inner  14  cells  of  the  cathode  endplate  inside  the  shank  were  designated  as  proton  'emis¬ 
sion*  points.  The  emission  at  those  points,  however,  was  completely  independent  of  the  adjacent  elec¬ 
tric  field  values.  A  consunt,  emitted  current  density.  J„  of  0.65  kA/cm^  was  imposed  for  the  1  kA 
beam  case  and  3.25  kA/cm^  for  the  5  kA  case  at  each  of  the  14  points.  Particles  of  both  species  were 
perfectly  absorbed  upon  hitting  any  of  the  boundary  surfaces,  including  the  surfaces  of  the  protruding 
cathode  shank. 

All  quantities  of  interest  in  the  device  were  monitored  by  the  computer  code's  extensive  diagnos¬ 
tics.  Electron  and  ion,  emitted  and  collected  current  densities  were  recorded  over  all  relevant  surfaces. 
Electrosutic  equipotential  contour  lines  and  magnetic  field  lines  throughout  the  device  were  plotted. 
Sample  values  of  all  three  components  of  the  magnetic  field  are  similarly  recorded.  Net  charges  and 
currents  were  periodically  listed.  Finally,  the  positions  of  sutistical  samplings  of  electrons  and  protons 
were  plotted  at  equilibrium. 

IV.  RESULTS 

First,  a  ’benchmark"  simulation  was  run  in  which  only  electfons  were  allowed  in  the  diode.  In 
agreement  with  theory  (see  Section  V),  a  net  electron  current  of  9.6  kA  was  mea«ure).  Of  that  toul 
amount.  2.8  kA  was  emitted  from  the  face  of  the  cathode  shank  tip  while  the  remaining  6.8  k.\  ori¬ 
ginated  from  the  outside  surface  of  the  shank.  Emission  from  the  inner  shank  surface  was  negligible. 
A  net  electron  charge  of  -1.86  x  10*  statcoulombs  was  found  to  be  in  the  system.  Figure  4  presents  a 
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plot  of  sample  electron  positions  at  equilibrium.  Of  note  is  the  waviness  of  the  outer  envelope  of  the 
electron  beam.  Similar  beam  behavior  was  observed  in  foilless  diode  simulations  conducted  by  R.  Jack- 
son*^  '^  and  others.  The  electrostatic  potential  contour  plot  for  the  device  modeled  in  simulation  is 
shown  in  Fig.  5.  The  beam  envelope  waviness  is  mirrored  in  these  equipotentials.  It  is  also  clear  from 
this  figure  that  the  electric  held  inside  the  hollow  cathode  is  negligible,  in  agreement  with  the  vanish¬ 
ingly  small  electron  emission  observed  there.  It  also  appears  that  electrons  leaving  the  cathode  are 
fairly  rapidly  accelerated  to  the  40%  equipotential  (0.2  MeV)  before  entering  a  virtual  "drift  region" 
about  3.S  cm  long  starting  about  2.0  cm  downstream  of  the  cathode  tip.  Electron  energy  and  velocity 
are  nearly  constant  in  this  region  as  would  be  the  case  in  an  actual  full-length  foilless  diode  drift  tube. 
It  is  therefore  that  region  where  beam  behavior  will  most  closely  mimic  a  realistic  device.  Finally,  the 
self-magnetic  field  lines  generated  by  the  e-beam  are  depicted  in  Fig.  6.  Imposed  over  this  self-held 
pattern  is  the  uniform  10  kG  background  field.  Up  to  about  Z  ~  8.5  cm,  the  beam  field  is  diamagnetic 
in  nature,  lowering  the  inside,  imposed  82  by  about  40-80  gauss  and  intensifying  the  held  near  the 
anode  tube  by  a  like  amount.  The  self-held  between  Z  8.5  and  10.0  cm,  where  electrons  receive  a 
hnal  0.3  MV  of  acceleration,  is  of  the  opposite  orientation.  However,  that  region  as  well  can  boast  of 
only  50-100  gauss  held  strengths. 

The  second  case  to  be  tested  was  that  involving  the  1  kA,  30  MeV  proton  beam  injected  in  the 
manner  described  in  the  previous  section.  The  net  electron  current  in  the  new  equilibrium  increased  by 
2  kA  to  1 1.6  kA.  Of  that  total  current,  4.3  kA  was  being  emitted  from  the  tip  (compared  to  2.8  kA  for 
the  electrons-only  case),  0.4  kA  were  coming  from  the  inner  shank  surface,  while  the  outer  shank  sur¬ 
face  emission  remained  practically  unchanged  at  6.9  kA.  1‘he  ion  beam  had  introduced  a  total  charge  of 
-1-4.06  X  lO’  statcoulombs  into  the  diode  while  the  electron  charge  had  increased  to  -2.11  x  10* 
statcoulombs.  Thus,  the  electron  beam  characteristics  of  the  diode  itself  had  changed  very  little.  In 
addition,  the  ion  beam  envelope  experienced  virtually  no  ballooning.  This  is  evident  from  Fig.  7  which 
presents  the  steady-state  ion  beam  envelope  superimposed  over  the  sample  electron  position  plot.  The 
qualitative  character  of  the  hollow  electron  beam  has  likewise  remained  essentially  unchanged.  The 
first  truly  mqjor  change  manifests  itself  in  the  equipotential  plot  of  Fig.  8.  Comparison  with  Fig.  S 
shows  the  distortions  caused  by  the  concentration  of  positive  space  charge  along  the  axis.  The  contours 
are  sucked  into  the  hollow  cathode  up  to  the  40%  potential.  This  accounts  for  the  tum-on  of  electron 
emission  there.  Also,  electrons  reaching  the  previously-mentioned  "drift  region"  between  Z  ~  4.5  to  Z 
—  8.0  cm  are  now  at  about  0.25  MeV  instead  of  the  0.20  MeV  of  the  electrons-only  case.  Finally,  Fig. 
9  shows  the  changes  made  to  the  self-magnetic  held  lines  generated  by  the  electron  beam.  The  magni¬ 
tudes  of  those  fields  in  the  drift  region  remain  unchanged  in  the  range  of  40  to  80  gauss,  but  the  orien¬ 
tation  of  the  field  has  reversed.  The  imposed  Bz  is  now  strengthened  inside  the  beam  and  weakened 
outside.  Although  this  effect  is  interesting,  the  small  field  strengths  make  it  unimportant  to  the  opera¬ 
tion  of  the  diode. 

For  the  third  and  final  test  case,  the  current  of  the  30  MeV  proton  beam  was  increased  to  5  kA. 
The  steady  state  electron  current  rose  to  19.7  kA  with  8.3  kA  from  the  shank  face,  4.7  kA  from  its 
inner  surface,  and  an  outer  shank  emission  current  almost  identical  to  the  previous  two  cases  at  6.7  kA. 
Of  course,  the  ion  charge  increased  by  the  same  factor  of  five  as  its  current  to  -4-2.04  x  10* 
statcoulombs.  In  partial  compensation,  the  net  electron  charge  had  grown  to  -3.32  x  10*  statcoulombs. 
Clearly  the  normal  operation  of  the  foilless  diode  has  been  seriously  disrupted.  Its  electron  current  has 
more  than  doubled.  Even  more  dramatic  is  the  destruction  of  the  ion  beam.  This  is  illustrated  by  the 
equilibrium  sample  particle  plot  of  Fig.  10.  The  proton  beam  envelope  has  expanded  halfway  through 
the  width  of  the  hollow  electron  beam  by  the  time  it  arrived  at  the  end  of  the  drift  tube.  Similarly  the 
electrosutic  potential  contours  shown  in  Fig.  1 1  no  longer  bear  much  resemblance  to  the  ion-free  case 
of  Fig.  5.  Electrons  now  experience  almost  80%  of  the  toul  diode  potential  increase  before  entering 
the  "drift  region."  There  no  longer  exists  any  effective  "potential  well"  through  which  the  proton  beam 
can  travel.  Finally,  the  self-magnetic  held  of  the  electrons  (see  Fig.  12)  can  no  longer  be  ignored.  Its 
drift  region  magnitudes  now  range  from  about  200  to  700  gauss  which  is  an  appreciable  fraction  of  the 
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imposed  10  kG  field.  As  in  the  1  kA  ion  beam  case,  this  self-held  strengthens  the  Bz  along  the  axis 
and  weakens  it  outside  the  e-beam  radius. 

The  major  results  of  the  above  three  simulations  are  summarized  in  Table  1.  The  statement  can 
be  made  that  the  injection  of  a  I  kA,  30  MeV  proton  beam  along  the  axis  of  the  foilless  diode  being 
studied,  does  not  seriously  disrupt  the  normal  operating  characteristics  of  the  device.  The  ii\jection  of  a 
5  kA,  30  MeV  ion  beam  does  cause  serious  disruption.  Correspondingly,  the  1  kA  beam  is  virtually 
uneffected  by  its  passage  through  the  diode,  while  the  S  kA  beam  undergoes  severe  ballooning.  In 
order  to  better  understand  these  macroscopic  results,  it  is  important  to  examine  the  modifications  made 
by  each  case  on  the  individual  current  density  profiles  of  each  case.  Figure  13  depicts  the  radial  profiles 
of  the  electron  current  density,  J„  emitted  from  the  tip  face  of  the  cathode  shank.  The  profiles  all 
have  a  characteristic  double-pronged  shape. They  show  that  for  this  diode  the  doubling  of  the  emitted 
currents  on  this  surface  is  predominantly  occurring  at  the  points  of  lower  radius.  Those  points  are 
closer  to  the  ion  beam,  whose  space-charge  has  grossly  distorted  the  electric  field  values  there.  From 
Fig.  14,  however,  it  is  obvious  that  the  electron  emission  from  both  the  inner  and  the  outer  shank  sur¬ 
faces  are  virtually  identical  for  the  electrons-only  and  the  1  kA  beam  case.  Since  some  70%  of  the  total 
diode  current  originates  on  the  outer  surface,  this  agrees  with  the  overall  close  similarity  between  the 
two  cases.  The  S  kA  beam  case  is  quite  a  different  matter.  Although  the  net  outside  shank  emission  is 
about  the  same  as  the  previous  cases,  the  axial  profile  of  J,  is  grossly  depressed  over  the  front  half  of 
the  shank  and  mildy  enhanced  over  the  rear  half.  This  is  a  manifestation  of  the  2  1 /2-fold  increase  in 
the  value  of  Bg  along  that  surface.  The  near-normal  amount  of  negative  charge  emitted  from  the  rear 
half  is  prevented  from  leaving  the  near  vicinity  of  the  shank  as  it  travels  toward  the  front  face.  This 
negative  space  charge  then  impedes  the  electron  emission  along  the  front  half  of  that  surface.  The 
most  dramatic  change,  however,  is  in  the  inner  shank  surface  emission.  For  the  5  kA  beam  case,  the 
emission  there  rivals  that  of  the  outside  surface.  For  the  other  two  cases  that  inner  surfaces  had  been 
effectively  non-emitting.  The  combined  emission  changes  are  all  reflected  in  the  radial  profiles  of  elec¬ 
tron  current  densities  collected  at  the  anode  face  (Z  ■■  10  cm).  These  are  shown  in  Fig.  15.  In  addi¬ 
tion  to  the  overall  enhancement  of  collected  current  from  case  to  case,  a  distinct  pinching  of  the  beam 
toward  the  central  axis  can  be  seen.  This  is  to  be  expected  due  to  the  increased  self-magnetic  field 
strength  but  the  magnitude  of  the  radial  displacement  of  the  mean  current  density  is  limited  by  ihe 
imposed  10  kG  axial  S-held. 

The  final  diagnostic  measured  the  radial  profiles  of  the  30  MeV  proton  beam  in  the  plane  located 
at  Z  —  7.5  cm,  well  inside  the  electron  beam  "drift  region."  Two  profiles  are  shown  for  the  1  kA  beam 
case  in  Fig.  16.  The  first  represents  the  initial,  imposed  beam  profile  a  mere  50  timesteps 
(Af  -  2  X  10"*^  sec)  after  the  beginning  of  the  simulation.  The  other  profile  shows  the  modest  spread 
of  the  beam  about  0.8  nano-seconds  later.  The  change  is  quite  negligible,  indicating  excellent 
confinement  and  stability  of  the  beam.  Once  again,  the  5  kA  case  is  dramatically  different.  The 
corresponding  profiles  for  that  high  current  beam  are  shown  in  Fig.  17.  Over  the  same  0.8  nanosecond 
period,  this  beam  has  expanded  from  its  initial  0.7  cm  radius  to  over  I.O  cm,  well  beyond  the  inner 
cathode  shank  radius  and  the  inner  radius  of  the  hollow  electron  beam.  That  configuration  is  clearly 
unstable. 

V.  THEORY  COMPARISON 

in  this  section  analytic  expressions  are  derived  and  compared  to  the  simulation  results.  The  two 
major  results  are  the  radial  ballooning  of  the  ion  beam  and  the  increase  in  the  electron  limit  current  as 
the  ion  current  is  increased.  Agreement  between  theory  and  computation  is  quite  good. 

The  radial  equation  of  motion  for  the  ions  decouples  and  only  the  ion  beam's  self  fields  determine 
the  ion  motion.  This  is  true  for  times  that  are  short  in  comparison  to  the  axial  gyro-period.  That  is. 

qB 

where,  fl.  -  — -  and  B.  is  the  axial  magnetic  field.  .Also,  for  j3^  <  <  1/2,  where  B  is 
y\f 
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the  ion  velocity  divided  by  the  speed  of  light,  then  changes  in  the  relativistic  mass  facto'  can  be 
neglected.  With  these  assumptions  the  equation  of  motion  for  the  ions  at  the  beam  surface  is  given  by 


qE,  -  qV^Bg 


(1) 


where;  q  and  M  are  the  ion  charge  and  mass,  E,  and  Bg  are  the  seif  electric  and  magnetic  fields  in  the 
ion  beam,  y  is  the  relativistic  mass  factor,  is  the  axial  ion  velocity,  and  r  is  the  radial  position  of  ions 
at  the  beam  surface.  Assuming  current  conservation  Eq.  (1)  reduces  to'^ 


(fir  k 
dfi  r 


(2) 


where;  k  —  - - T'  ^  the  ion  beam  current,  and  en  is  the  dielectric  constant  of  free  space. 

2ireoA/Kjy^ 

Equation  (2)  can  be  solved  in  series  form‘*  to  give. 


,3, 

*5  (2/i  +  l)«!  /-o 

where;  tq  is  the  initial  beam  radius  at  r  —  0  and  r  is  the  final  radius  at  time  t.  In  Table  2,  a  comparison 
is  made  between  the  code  time,  r,  to  reach  the  final  radius,  r,  and  the  time  calculated  by  Eq.  (3).  This 
comparison  is  made  for  30  MeV  protons. 


It  is  clear  from  the  comparison  that  the  self  fields  of  the  ion  beam  determine  its  motion.  Of  the 
self  fields,  the  electric  field  dominates  in  the  expansion.  This  result  does  not  limit  the  maximum  ion 
current.  Since,  the  ion  beam  expands  to  the  diameter  of  the  electron  beam,  it’s  electric  field  is  reduced 
thus  allowing  for  more  ion  current  to  flow.  Increasing  the  ion  current  then  increases  the  electron  limit 
current  which  in  turn  allows  for  an  increase  in  the  ion  current,  etc. 


In  addition  to  the  problem  of  ion  beam  ballooning,  the  simulation  also  addressed  the  major  ques¬ 
tion  of  changes  in  the  diode's  electron  current  due  to  the  presence  of  the  ion  beam.  The  maximum 
electron  current  is  determined  by  the  potential  disuibution  along  the  electron-emitting  cathode  surfaces. 
When  ions  are  present  the  potential  is  raised  and  hence  the  limit  current  is  increased. 

In  Fig.  18  electrons  and  ions  are  moving  in  the  same  axial  direction.  A  large  external  magnetic 
field  is  applied  in  the  direction  of  particle  motion,  z.  This  field  is  large  enough  to  magnetize  the  elec¬ 
trons  but  not  the  ions.  Furthermore,  the  assumptions  are  made  that  (a)  J,  and  are  uniform,  (b) 
\Bg,\  and  «  Ifl,!,  (c)  <  Ifljlc,  and  (d)  Electron  velocity  shear  can  be  neglected  (where  7* 

and  /  are  the  electron  and  ion  current  density;  Bg/,  Bg^  and  B^  are  the  ion,  electron  and  external,  mag¬ 
netic  fields;  Ert  is  the  radial  ion  electric  field;  c  is  the  vacuum  speed  of  light) . 


From  Gauss'  law  the  potential  inside  the  electron  beam,  region  3,  (see  Fig.  20)  is  given  by  the 
expression,'’ 


<i(r) 


Ij 

2ireoK 


In  (/?/r) 


4ireoFj 


’’i  ~ 

- j-  2  In  R/ri  - 

rj  -  rf 


In  (r^/r) 


(4) 


where  /;  and  /,  are  the  ion  and  electron  current,  FJ  and  are  the  ion  and  electron  velocity,  and  eo  is 
the  permitivity  of  the  free  space.  The  quantity,  <^,  may  be  eliminated  from  Eq.  (4)  by  using  conserva¬ 
tion  of  energy, 
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mc^yo  is  the  initial  electron  energy,  m  and  e  are  the  electronic  mass  and  charge,  and  /3,  is  the  ion  velo¬ 
city  normalized  to  the  vacuum  speed  of  light. 


By  inspecting  the  first  term  in  Eq.  (8)  it  is  clear  that  the  presence  of  the  ion  current  raises  the 
final  electron  energy.  This  energy  increase  plays  only  a  small  role  in  the  limit  current  expression.  The 
biggest  ion  effect  arises  in  Eq.  (4)  where  the  potential  is  directly  increased  by  an  increase  in  the  ion 
current. 


The  magnitude  of  the  limiting  electron  current  of  Eq.  (6)  varies  with  the  radial  position,  r.  In 
general  Eq.  (6)  should  be  integrated  from  to  rj  to  give  a  current  that  is  weighted  by  the  changing 
potential  across  the  beam.  However,  since  the  equation  itself  is  only  a  simple  approximation,  seeking 
such  extra  "precision"  would  be  pointless.  Instead,  Eq.  (6)  will  be  evaluated  at  different  radii  and  the 
effects  discussed.  An  upper  bound  to  the  limit  current  can  be  obtained  by  setting  r  ^  then  Eq.  (6) 
becomes 


^<21^  I-  (1  - 

l\nRlri  fit 


This  gives  an  upper  bound  to  the  limit  current  since  all  particles  are  experiencing  the  maximum  poten¬ 
tial  near  the  wall.  A  comparison  is  made  in  Table  3.  If  Eq.  (6)  is  evaluated  at  r  —  rj  then  the  result  is. 


Note  that  when  no  ions  are  present,  this  will  give  a  lower  bound  to  the  limit  current  since  ail  of  the 
particles  are  exposed  to  a  potential  minimum.  Table  4  shows  a  comparison  of  Eq.  (10)  with  the  code 
results.  As  expected,  the  theoretical  prediction  is  well  below  the  code  findings  for  the  case  with  no  ion 
beam  present.  For  the  two  cases  with  ion  injection,  however,  the  addition  of  positive  space  charge 
inside  r  —  /-j,  elevates  the  potential  there  to  such  a  degree  that  for  the  S  kA  case  it  actually  exceeds 
that  observed  in  Table  3  for  /■  -  ri-  It  is  clear  from  both  Table  3  and  Table  4  that  the  comparison 
between  simulation  and  theory  is  quite  good,  considering  the  theoretical  approximations  that  were 
made.  Thus  a  scaling  law  for  the  electron  limit  current  dependence  on  the  ion  current  can  be  inferred. 
This  effect  can  have  a  major  impact  on  wave  phase  velocity  control  for  the  space  charge  wave  accelera¬ 
tor.  For  example  if  the  injection  current  is  fixed  a  change  in  the  ion  current  will  shift  the  wave  phase 
velocity  thus  removing  ions  from  wave  synchronization.  However,  at  high  yo  and  for  modest  ion 
currents  the  effect  is  small.  For  example,  if  /,//,  -  0.1  and  yo  -  7.0, 
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Thus,  the  two  major  code  results  have  been  explained.  Analytic  expressions  have  been  derived 
for  predicting  these  results  and  good  agreement  was  found  between  theory  and  code  results.  The 
mechanism  for  ion  ballooning  has  been  shown  to  be  due  primarily  to  the  self  electric  field  of  the  ion 
beam.  Increases  in  the  electron  limit  current  are  shown  to  be  associated  with  an  increase  in  the  ion 
current.  As  discussed,  a  coupling  between  the  ion  and  electron  currents  can  negatively  impact  on  ion- 
wave  synchronization  for  the  space  charge  wave  accelerator. 

VI.  CONCLUSIONS 

The  numerical  simulations  presented  here  give  the  first  documented  evidence  that  a  high  current, 
high  energy  light  ion  beam  can  be  injected  down  the  central  axis  of  a  conventional  foilless  diode 
without  either  seriously  disturbing  the  operating  characteristics  of  the  diode  or  causing  significant  disr¬ 
uption  of  the  ion  beam  itself.  For  the  4-5  GW  foilless  diode  studied  and  for  a  proton  energy  of  30 
MeV  a  beam  current  of  1  kA  can  easily  be  tolerated.  Increasing  the  ion  beam  current  to  S  kA  leads  to 
a  general  break-down  of  the  beam-diode  system.  The  precise  current  value  between  1  kA  and  5  kA  at 
which  system  stability  ceases  remains  to  be  found.  At  y >7  a  ten  percent  ion  current  will  cause  less 
than  a  one  percent  variation  in  beam  energy.  This  appears  acceptable  for  collective  wave  accelerators. 
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Table  1  —  Diode  Currents  and  Charges 


Case  1 

Case  2 

Case  3 

/ion  (kA) 

0 

1 

5 

^eieciron 

9.6 

11.6 

19.7 

A,p  ^kA) 

2.8 

4.3 

8.3 

^outside 

6.8 

6.9 

6.7 

Anside  (kA) 

— 

0.4 

0.7 

Qe 

(xlOOO  statcoul) 

-18.6 

-21.1 

-33.2 

a 

(xlOOO  statcoul) 

0 

+4.06 

+20.4 

Table  2  —  Comparison 
of  Code  and  Analytic  Calculated  Times 
to  Reach  Final  Radius 


/ion  (kA) 

1 

5 

rg  (cm)  initial 

.725 

.725 

r  (cm)  final 

.83 

1.07 

code  time,  tins) 

.8 

.8 

time  from  Eq.  (3),  r(«s) 

■BO 

%  Relative  deviation  in 
times 

H 

11 

Table  3  —  Comparison  of 
Code  and  Analytic  Limit  Currents. 
Potential  evaluated  at  r  ■«  ri. 


/ion  (kA) 

0 

1 

5 

/electron  ^kA)  from 

Eq.  (9) 

11.03 

.  .. 

13.49 

24.13 

^d€ciron  ^^A)  from 

Code 

9.6 

11.6 

19.7 

%  Relative  Deviation 
from  Eq.  (9) 

12.7 

Table  4  —  Comparison  of 
Code  and  .Analytic  Limit  Currents. 
Potential  evaluated  at  r  -  rj , 


/ion  (kA) 

0 

»■  .1- 

^electron  ^kA)  ffom 

Eq.  (10) 

m 

11.01 

27.85  { 

1 

/el«:.ron  <kA)  from 

Code 

9  6 

11.6 

19.7  j 

.  i 

%  Relative  Deviation 
from  Eq.  HO) 

31.5 

29 
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M4ve  Acceleration  Wave  Growth 


V  =  0.5MV 
I  =  1.4kA 
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Fig.  1  —  The  converging  guide  accelerator 


Fig.  2  —  The  experimental  diode  geometry 
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Fig.  3  —  The  numerical  diode  geometry 
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Fig.  4  —  Sample  electron  positions  in  equilibrium  diode 
with  no  ions  present 
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Fig.  5  —  Mr.  z)  for  electrons— only  diode 
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Fig,  6  —  Self-magnetic  field  for  electron— only  diode 


diode  with  1  kA. 
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Fig.  15  —  /,(r)  collected  at  anode  endplate 


Fig.  16  —  Profiles  of  /,(rl  for  1  kA, 
30  MeV  proton  beam  case 
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Fig.  17  —  Profiles  of  J,{r)  for  I  kA. 
30  .MeV  proton  beam  case 
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